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Abstract. We present numerical simulations of three-dimensional compressible magnetoconvection in a rotating 
rectangular box that represents a section of the solar convection zone. The box contains a convectively unstable 
layer, surrounded by stably stratified layers with overshooting convection. The magnetic Reynolds number, Rm, 
is chosen subcritical, thus excluding spontaneous growth of the magnetic field through dynamo action, and the 
magnetic energy is maintained by introducing a constant magnetic field into the box, once convection has attained 
a statistically stationary state. Under the infiuence of the Coriolis force, the advection of the magnetic field results 
in a non- vanishing contribution to the mean electric field, given by {u x b). From this electric field, we calculate 
the a-effect, separately for the stably and the unstably stratified layers, by averaging over time and over suitably 
defined volumes. From the variation of a we derive an error estimate, and the dependence of a on rotation and 
magnetic field strength is studied. Evidence is found for rotational quenching of the vertical a effect, and for a 
monotonic increase of the horizontal a effect with increasing rotation. For Rm ~ 30, our results for both vertical 
and horizontal a effect are consistent with magnetic quenching by a factor [1 + Rm{Bo/ Bcq)^]~^ ■ The signs of 
the small-scale current helicity and of the vertical component of a are found to be opposite to those for isotropic 
turbulence. 
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1. Introduction 

Magnetic fields are observed on a wide variety of cosmi- 
cal bodies, among which are planets, stars and galaxies. 
With the exception of a few types of objects whose mag- 
netic fields are thought to be frozen-in relic fields, cosmical 
magnetic fields are attributed to dynamo action. Dynamo 
theory concerns the generation of magnetic fields in elec- 
trically conducting fluids. In stars and planets, dynamo 
action is the result of an interplay between convection 
and rotation. The present simulations focus on dynamo 
action in late-type stars, which are characterized by an 
outer convective zone on top of a stably stratified radia- 
tive interior. Many late-type stars are magnetically active, 
and some exhibit magnetic cycles, as does the Sun. 

A successful solar model was first proposed by Parker 
(1955) who recognized that shear and hehcity cooper- 
ate in such a way that the mean magnetic field oscil- 
lates in a migratory manner; later Yoshimura (1975) has 
shown that the field generally propagates along the sur- 
faces of constant angular velocity. Since then, magnetic 
cycles and butterfly diagrams have been produced by 
mean-field models using spherical geometry (Steenbeck & 



Krause 1969), with depth-dependent magnetic diffusivity 
(Roberts & Stix 1972), and with the character of an in- 
terface wave at the base of the convection zone (Parker 
1993, Charbonneau & MacGregor 1997). As an alternative 
to the helical effect of convection (the a effect), the sys- 
tematic tilt of flux tubes that become unstable and erupt 
to the solar surface has been employed as a regenerative 
agent for the mean poloidal field (Leighton 1969; Durney 
1995, 1997; Schmitt et al. 1996; Dikpati & Charbonneau 
1999), with equal success regarding the existence of Sun- 
like cycles and butterfly diagrams. In all cases, however, 
the regenerative terms in the mean-field equations were 
based either on approximations such as first-order smooth- 
ing, or on the limited available knowledge about the be- 
havior of magnetic flux tubes in the convection zone; the 
review of Stix (2001) summarizes some of the current 
problems. 

In principle, the problem can be attacked by direct nu- 
merical integration of the equations of magnetohydrody- 
namics (MHD) for the st ar's co nvection zone (for a review, 
see e.g. Nordlund et al. 1994). But severe problems are 



encountered. Firstly, a stellar convection zone is a highly 
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turbulent plasma (Re > 10^^), so that the kinetic energy 
spectrum encompasses a very wide range of length scales. 
If these were all to be resolved, a prohibitive number of 
mesh points would be required. Secondly, since the mag- 
netic Reynolds number is also large (Rm > 10*), a similar 
argument holds for the length scales of the magnetic field. 
For these reasons three-dimensional MHD simulations of- 
ten have been restricted to a rectangular box that rep- 
resents only a small section of a stellar convection zone. 
In the present work we do this, too. Thirdly, the Prandtl 
number, the ratio of the kinematic viscosity to the radia- 
tive diffusivity, is very small in the solar convection zone 
(Pr « 10~^ — 10~^), so that the flow can vary on much 
smaller scales than the temperature. This is known to have 
consequences for the topology and te mpor al behavior of 
conve ctive turbulence (Cattaneo et al. 1991 ; Brandenburg 
et al. 1996| ; Brummell et al. 199(j ), but in the present sim- 
ulations, we ignore such effects, and set Pr = 1. Fourthly, 
the Mach number. Ma = u/cg, decreases with depth in the 
solar convection zone from about 0.1 in the photosphere 
to less than 0.001 near the base. The corresponding high 
sound speed in the lower convection zone necessitates the 
use of very small time steps in order to fulfill the CFL 
condition for compressible hydrodynamics {At < Ax/cg), 
while it is often desirable to continue simulations for at 
least several convective turnover times. One possibility to 
circumvent this problem is to adopt the anelastic approx- 
imation whereby sound waves are excluded, so that larger 
time steps can be used (Ogura & Charney 1962; Lantz & 



Fan 1999). Here we maintain full compressibility, and at- 
tain a mean Mach number of the order 0.1; in a subsequent 
paper we shall study effects of a smaller Mach number. 



Brandenburg et al. ( |1990[ ) extended a 3D hydrodynam- 
ical code to the case of magnetoconvection, including the 
effect of rotation. Subsequently, a spontaneous dynamo 
inst abilit y was observed in th e simu lations (Nordlund ct 
al. 



1992 



1996). Recent simulations 



Brandenburg et al. 
of isotropically forced helical turbulence (Brandenburg 
2001 ) have verified the existence of large scale dynamo 
action, and it was possible to identify this as the re- 
sult of an a effect (in the sense of a non-local inverse 
cascade). However, the large scale field generated pos- 
sesses magnetic helicity and, since for closed or periodic 
boundaries the magnetic helicity can change only resis- 
tively, the growth of the large scale field is slowed down 
as the magnetic Reynolds number increases. This trans- 
lates inevitably to a magnetic-Reynolds-number depen- 
dent a effect and turbulent magnetic diffusivity, as sug- 
gested by Vainshtein & Cattaneo (1992). The argument of 
Vainshtein & Cattaneo is however only phenomenological, 
not based upon the fundamental concept of magnetic he- 
licity conservation, and hence their conclusion that strong 
large-scale fields are impossi ble is not borne out by the 
simulations of Brandenburg (2001). Furthermore, in that 
paper it was shown that the a obtained by imposing a 
magnetic field is indeed a reasonable approximation to 
the a that results naturally even if no field is imposed. 



In the present simulations we do not calculate box dy- 
namos, but rather concentrate on the dynamo coefficients 
that occur in mean-field theory. Therefore, all simulations 
are performed with an imposed magnetic field, and the 
magnetic Reynolds number is chosen subcritical, so that 
the magnetic field is not self-sustained. We think that 
progress in stellar dynamo theory can be made through 
a combination of exact MHD simulation and mean-field 
dynamo theory. Furthermore, the existence of systematic, 
large-scale magnetic fields and stellar cycles with peri- 
ods far in excess of convective time scales suggests that 
a mean-field description is possible in some form. Also, we 
argue that mean-field theory (in the wide sense that in- 
cludes the above-mentioned models such as Leighton's) is 
the only current model that reproduces large-scale mag- 
netic fields and cycles, including such an outstanding fea- 
ture as the solar butterfly diagram, even though the di- 
verse approximations made for calculating the dynamo co- 
efficients may not be valid under stellar conditions. The 
case of self-excited dynamo action would of course be 
very interesting too, especially in view of the question of 
whether a-quenching depends on the magnetic Reynolds 
number, as suggested by Vainshtein & Cattaneo (1992). 
However, we postpone this until a later study, partly be- 
cause the measurement of the dynamo coefficients in the 
presence of a large-scale field that is different from the 
imposed one is less straightforward. 

The present investigation is limited to the regime 
where the coefficients are determined directly by the flow, 
and we do not address the question of whether the large- 
scale magnetic field of late-type stars is in fact generated 
by such an ordinary a effect in the convection zone, or 
by a magnetic instability in an underlyin g stab ly strat- 
ified layer (e.g., Brandenburg & Schmitt 1998 ). We do 
include a stably stratified layer with overshooting convec- 
tion, but its main purpose here is to provide realistic con- 
ditions for the flow in the unstable region. In most runs, 
the strength of the imposed field is set to a value which 
amounts to typically 2% of the equipartition field with re- 
spect to the kinetic energy density. During the subsequent 
evolution the field strength also remains small compared 
to the equipartition value. Furthermore we shall explore 
the influence of the imposed magnetic field by increas- 
ing its strength up to values somewhat in excess of the 
equipartition value. Therefore the results may have rele- 
vance for the solar convection zone, where the magnetic 
field is no stronger than the equipartition value. The a ef- 
fect produced in the galactic gas by supernova explosions 
has been calculated by Ziegler et al. (1996), in a similar 
spirit as in the present work. 

The relevance of the transport coefficients for mean- 
field dynamo theory becomes clear from the equation for 
the mean magnetic field. 



d{B) 
dt 



= Vx [(J7)x(B) + (mx5) -tjS/x{B) 



(1) 



where (B) is the mean magnetic field, (U) is the mean 
fiow, u = U — (U) is the fiuctuating component of the 
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flow, and {u x b) is a contribution to the mean elec- 
tric field, sometimes referred to as the electromotive force. 
The dynamo coefficients appear in an expansion of this 
mean electric field in terms of spatial derivatives of the 
mean magnetic field. In general they can be represented 
as kernels of an integral equation (Brandenburg & Sokoloff 



2001). In the simplest case, the coefficients are treated as 
local tensors, which leads to 



{u X b), aij{Bj) + f3^jk^j{Bk) + 



(2) 



The first term in the expansion is the a effect. The alpha 
tensor is a pseudo tensor that exists only in non mirror- 
symmetric flows, as they occur in stellar convection. For 
the solar dynamo, the main significance of the a effect lies 
in generating the poloidal mean magnetic field, which is 
achieved predominantly by a^^. Toroidal fields are gener- 
ated by the strong differential rotation that exists in the 
tachocline near the base of the solar convection zone. In 
other solar-type stars and fully convective stars though, 
differential rotation could be weak (Kiiker & Stix 2001), 
so that the dynamo becomes of the a^-type. In that case, 
other components of the alpha tensor, such as ct^|. , are 
important for maintaining the toroidal magnetic field. 

The remainder of the paper is structured as follows. 
After introducing the model and the equations, we focus 
on the dependence of a on rotation. In the following sec- 
tion, effects of the strength and orientation of the imposed 
field are studied. In the final section, a discussion of the 
results is presented. 

2. The Model 

The model used in the simulations is the same as that of 
Brandenburg et al. (1996). The simulation domain con- 
sists of a rectangular box which is defined on a Cartesian 
grid, with x and y denoting the two horizontal coordinates 
(corresponding to latitude and longitude in spherical co- 
ordinates), and z denoting depth, respectively (Fig. 1). In 
most cases, convectively stable layers are included below 
and above the unstable layer. The upper layer (region 1) is 
stabilized by a cooling term in the energy equation, which 
leads to an almost isothermal, highly stable stratification. 
The lower stable layer (region 3) represents an overshoot 
zone, whose thickness is chosen such that overshooting 
plumes do not reach the lower impenetrable boundary. 
The box rotates about an axis fl, that makes an angle 
with the z-axis. Hence 9 = corresponds to a situation 
where the box is located at the south pole of the Sun. In 
the present paper we treat only this special case. 

The governing equations are those describing magnetic 
induction, mass continuity, and the balance of momentum 
and energy: 



dB 

'dt 

dp 

'dt 



Vx{UxB)+rjV^B , 



(3) 
(4) 




Fig. 1. Geometry of the simulation domain. The layer 
boundaries are at zi = —0.15, Z2 = 0, Z3 = 1, and 
Z4 = 2.85, except for run 10 that has no stable regions, so 



that zi 
is 4. 



dU 

de 



Z2 = and Z3 



Z4 



1. The horizontal extent 



-\7p + pg + JxB-2pnxU + 2iy\/-pS, (5) 
-pV-U + \'-KVe + 2iypS'^ + HarjJ^ + Q , (6) 



where «/ = V x B//iq is the current density, 77 is the mag- 
netic diffusivity, v is the kinematic viscosity, 7 = Cp/Cy 
is the adiabatic index, and k is the radiative conductivity. 
In all simulations, rj and 1^ are taken constant, and k is a 
prescribed function of depth. The stress tensor S is given 

by 



1 / dUi dU, 



1 



2 I dxj dxi I 3 



(7) 



and stands for J^ij ^Ij ■ The equation of state is that of 
an ideal gas, i.e.. 



The term Q, given by 
Q = -<^Qf{z) (e - ei) , 



(8) 



(9) 



represents cooling or heating, depending on whether the 
internal energy density exceeds ei or falls below it, respec- 
tively; ctq represents the cooling/heating rate. The depth- 
dependent function / ensures that Q is non- vanishing only 
in region 1, and /(Z2) = 0. The inclusion of this term leads 
to the presence of a highly stable, thin overshoot layer, 
thereby providing a realistic upper boundary condition for 
the convection zone. In the horizontal directions, periodic 
boundary conditions are imposed. The upper and lower 
boundaries of the domain are impenetrable and stress-free, 
and the horizontal components of the magnetic field vari- 
ation are set to zero. At the lower boundary, the energy 
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flux is prescribed, and at the top boundary the internal 
energy (which is proportional to the temperature) is fixed. 









dux duy 








dz dz 


z = 


Zl, Z4 


< 


Bx — Bqx , 








. By = Boy , 


z = 


Zl 


e 


= ei , 






de ( 


z = 


Z4 


dz \dz / A 



(10) 



In the actual runs we set Boy = 0; Bqx is cither or the 
imposed horizontal field. 

All quantities are made dimensionless by setting 



d = po = g = m = = 1 



(11) 



where po is the initial density at a depth z — z^, the 
bottom of the unstable layer. Thus length is expressed in 
terms of d = z^ — Z2, the thickness of the convectively 
unstable zone (region 2). It follows that time is measured 
in terms of ^ dj velocity in terms of \fgd^ the magnetic 
field strength in terms of ^/popogd, and entropy in terms 
of Cp. 

On several occasions we shall consider the energy bal- 
ance, which is governed by a conservation law, 



dt 



(12) 



where etot — e + |itp/2 + \B\'^/2poP is the total specific 
energy, and the fiuxes, Fi, are given by 



-t'conv = jpeu, 
Fki„ = p\u\^u/2, 

-Frad = -KVe, 

-Fvisc = -2vpS : u, 
Fcm = ExB/po- 



(13) 
(14) 
(15) 
(16) 
(17) 



Apart from the geometry, boundary conditions, and ini- 
tial conditions, solutions are governed by the following 10 
independent dimensionless parameters. The Prandtl num- 
ber and the magnetic Prandtl number are defined as 



Pm ~ v/rj , 



(18) 



where xo — 1^2/1 Pa denotes a reference value of the ther- 
mometric (radiative) diffusivity for region 2. The param- 
eter ^0 determines the pressure scale height at the top of 
the box, ^0 = Hp{zi)/d: 



Co = (7- ^)ei/gd. 



(19) 



The initial thermal structure of the (maximally) three re- 
gions is characterized by the radiative temperature gradi- 
ents. 



'dlnT 



. d\np ' i 



(« = 1,2,3) 



(20) 



In addition, one often employs the polytropic index, m = 
(1 — V)/V. The adiabatic temperature gradient is given 
by Vad = (7 ^ l)/7- A measure for instability is provided 
by the superadiabaticity, 5i — Vad, which is positive 
in an unstably stratified medium. The Rayleigh number 
is defined as 



Ra = 



I'XoHph ' 



(21) 



where Hp^ = S,od + 0.5c?/(m2 -I- 1) is the pressure scale 
height in the middle of the unstable layer, as can be shown 
using the hydrostatic equilibrium (p5|). The parameters 82 
and Hph refer to the unperturbed stratification, i.e. that 
before the onset of convection. The Rayleigh number is a 
measure of the strength of convection compared to that 
of viscous and thermal (radiative) dissipation, as can be 
seen by writing Ra = ivisc^rad/^conv, where ^visc = rf^/i^, 
trad = d'^/xo, and t^^^^ = Hph/g62. Alternatively, one 
may express Ra in terms of the entropy gradient, ds/dz = 
Cp 8 1 Bp. According to the Schwarzschild criterion, a posi- 
tive value for (5, i.e. Ra > 0, signifies instability. In reality, 
Ra must exceed a finite threshold value for convection 
to set in. It should be noted that the local value of the 
Rayleigh number, d'^gS / (vxHp), varies with depth within 
the unstable layer, because x. Hp and, as a result of con- 
vection also 5, are z-dependent. Typically, the local value 
is several times smaller than Ra, mainly because 5 can be 
strongly reduced by convection. The Taylor number. 



Ta = {2nd^/vf 



(22) 



measures the importance of rotation relative to viscous 
dissipation. Finally, do represents the rate at which inter- 
nal energy is lost from the upper stable layer, and 6 is the 
angle between the rotation vector and the z-axis. 

All other parameters are secondary. The Coriolis num- 
ber, or inverse Rossby number, measures the importance 
of the Coriolis force and is defined as 



Co = 2nT , 



(23) 



where r — ^/urms is the turnover time. The correlation 
length, is taken to be d in the unstable region. The 
Chandrasekhar number. 



2 j2 



Ch 



mB^d 
ATrpoiyi] ' 



(24) 



measures the strength of the imposed magnetic field. In 
the initial state, the z-component of the radiative en- 
ergy flux, Fiad,2 = —ude/dz, is assumed to be constant 
throughout the domain. This determines the radiative 
conductivities in the three regions according to Ki/K2 = 
(rrii -|- l)/(m2 -l- 1). In fact, n is turned into a smooth 
function of depth by allowing it to change continuously 
across thin intermediate layers between the three regions. 
An approximate initial stratification, unperturbed by con- 
vection, is calculated on the assumption of hydrostatic 
equilibrium, and this is done iteratively until the condi- 
tion p{z = Z3) = Po is satisfied. It is also assumed that 
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region 1 is cooled efficiently enough to become isothermal. 
The result is then basically a smoothed version of 



e(zi) 
6(23) 

P(2l) 



Z — Z\ 


(mi 


+ l)i?„(zi)- 






(m2 






Z- Z'i 


(ms 


+ l)i?p(^3)J 




Z — Z\ 


(mi 


+ l)Hp{zi). 




Z - Z2 


(m2 


+ l)Hpiz2). 




z — z^ 



(ma + l)i?p(z3) 



(i) , 

(ii) , 

(iii) , 

(i) , 

(ii) , 

(iii) , 



(25) 



(26) 



where (i) ■ ■ ■ (iii) stand for the three regions, zi < z < Z2, 
22 < 2 < Z3, and z^ < z < Zi, respectively. The ac- 
tual initial stratification in region 2 is calculated nu- 
merically using the mixing-length formalism of convec- 
tion. The advantage of this approach is that it reduces 
the amount of time required for relaxation to a fully 
convective state. The internal energy density at the top 
equals ei = ■CoWadffrf- Using (25), it is easily shown 
that {de/dz)i = rui^g / (m^ + !)■ The radiative, kinetic 
and magnetic diffusivities follow from Eqs. (^8|) and (^l|). 
Reynolds numbers are defined as 



Re 



Rm ^ u,.^sd/ri, 



(27) 



where ttrms is the rms velocity defined in a suitable way 
(e.g., by averaging over time and over a partial volume of 
the box). 

We employ a finite difference scheme, according to 
which spatial derivatives are calculated with 6th-order ac- 



curacy (Lele 1992 ). Time-stepping is done using a third- 
order Hyman predictor-corrector method. Table |l| gives 
a list of the parameters used for a first series of runs in 
which the influence of rotation is investigated by varying 
the Taylor number. 

3. The a effect 

Essentially, the a effect amounts to a proportionality be- 
tween the mean electromotive force, (uxb), and the mean 
magnetic field. Magnetoconvection in a rotating, stratified 
medium provides the necessary anisotropics for a non-zero 
a. In the present case, with the rotation axis chosen par- 
allel to the direction of gravity, the most elementary, non- 
isotropic form of the expression for the a effect is one that 
distinguishes vertical and horizontal components: 



{uxb) = ay{Bz) + aii{Bi 



(28) 



The components an — a^x — ctyy and ay — ctzz cor- 
respond to a^^ and arr in spherical coordinates, respec- 
tively. We recall that a^^ plays the dominant role in the 
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Fig. 2. Averaging procedure (run 6). Top: at each time 
step, ay is calculated as a function of depth by averaging 
over the horizontal coordinates (grey curves). The thick 
drawn curve is the time average, and the surrounding thin 
drawn curves indicate the error in the mean. Bottom: ay 
is calculated separately for two partial volumes, namely 
Q!v,unst for the depth range —0.15 < z < 1, and ay.stab for 
the range 1 < z < 1.5. The error in the time averages is 
determined on the basis of the number of coherence times 
covered by the simulation. All calculations are for the solar 
south pole. 



generation of the poloidal magnetic field. The two com- 
ponents of a are extracted from the simulations by con- 
structing ay = (uxb) ■ (Bz) /\{Bz)\'^ and an = (uxb)- 
(Bh)/| (-Bh) In practice, (uxb) in these expressions is 
replaced by {UxB), because the difference, (U) x {B), 
was found to be negligible in every case. Depending on 
which component of a, i.e. ay or an, is to be determined, 
a magnetic field is imposed in the z-direction or in the 
cc-direction, respectively. 

3.1. Statistics 

It is well known from other numerical simulations that 
a is an extremely noisy quantity, so that one should do 
as much averaging as possible in order to optimize the 
statistics. 

For this paper we have chosen two averaging proce- 
dures. The first, which is considered in the present sec- 
tion, consists of an average over the horizontal coordinates 
and over time (Fig. |2[ upper panel). The time average 
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Table 1. Run parameters. For the geometry, see Fig. 1. The imposed magnetic field is Bq = BqEz, or BqBx, depending 
on whether ay or qh is determined. In all runs, the polytropic indices are given by mi = 0, TO2 = 1 and 777,3 = 3, 
the mesh size is 50'^ (except for run 10, where it is 30"^); Bq = 10^"^, Ra = 5 • 10^, Pr = Pm = 1 (i.e. Rm = Re), and 
^0 = 0.2. It follows that ly ~ 0.00211, and Ch — 0.018. The parameters Ma, Re and Co are averages over time and 
over the unstable layer; T is the time interval over which the averaging is performed. 



run 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Ta 


2 


10 


30 


100 


300 


2000 


10000 


20000 


50000 


10000 


T (ay) 


206 


822 


206 


841 


698 


545 


940 


362 


376 


554 


T (oh) 








564 


564 


456 


508 


346 


290 


521 


Ma 


0.089 


0.088 


0.088 


0.084 


0.091 


0.089 


0.080 


0.070 


0.056 


0.10 


Re 


34 


34 


34 


34 


35 


34 


30 


30 


20 


38 


Co 


0.042 


0.093 


0.16 


0.30 


0.50 


1.3 


3.3 


5.4 


11 


2.6 



is initiated from the moment when the magnetoconvec- 
tion has attained a statistically stationary state, i.e. when 
the kinetic and magnetic energy densities are more or less 
constant. The second is an additional average over two 
suitably chosen depth ranges (Fig. ||, lower panel): since 
the a effect depends on depth, and is expected to change 
sign near the bottom of the unstable layer, the averaging 
over depth is performed separately for the unstably and 
the stably stratified layers. This second procedure will be 
used mainly in later sections. 

The statistical deviation of a from its mean value gen- 
erally is large, often larger than the difference between 
that mean and zero, as illustrated in Fig. ||. Therefore, 
extended simulations are necessary in order to obtain a 
significant result. In the present paper the error estimate 
is based on the assumption that a single simulation can be 
divided into a number of time intervals containing inde- 
pendent realizations of the small-scale flow and magnetic 
field. The length r of those time intervals should then ex- 
ceed the correlation time, which is ~ 20 (in dimensionless 
units), as determined from the width of the autocorrela- 
tion function. To be on the safe side, we smoothed curves 
of Q!v,unst(i) and Q!v,stab(i) for several runs by performing 
a box average, and then calculated the autocorrelation 
function. This procedure yields a typical coherence time 
T = 50, which is several times longer than the value de- 
rived from the unsmoothed curves. 

Our procedure of estimating errors is more primitive 
than that applied earlier to results of numerical simula- 
tions (Pulkkinen et al. 1993). Moreover, the errors of the 
mean values determined in this manner can only be crude 
estimates because firstly we cannot be sure that we deal 
with Gaussian statistics (although inspection of the exam- 
ples of Fig. H suggests that a varies in a nearly symmet- 
ric manner around its mean), and secondly the number 
T/t of time intervals is not very large. Nevertheless we 
think that the procedure indicates whether or not the ob- 
tained mean values are significantly different form zero. 
Test simulations with identical parameters but different 
initializations confirm this indication. Of course, the re- 
sulting error bars do not imply anything about the phys- 
ical assumptions made. In particular, we cannot expect 



that models based on different assumptions, such as ours 
and that of Riidiger & Kitchatinov ( 1993 ), will yield re- 
sults that lie within error bars obtained from statistics 
(these bars shrink to zero for T 00). In contrast, we 
must see whether we can find qualitatively similar or dis- 
similar results from different models, especially as none of 
the models yet meets the real Sun. 

3.2. Relations between a and turbulence properties 

The existence of the a effect is attributed to the helical 
nature of convective flows in a rotating medium. Figure 
^ shows ay, an and several quantities that, as is shown 
below, have a relation to the a effect. Three different runs 
are selected in order to demonstrate the effects of rotation 
and boundary conditions. 

The a coefficients are depth-dependent, and undergo 
a sign change within the unstable layer. For Co ^ 4, the 
vertical alpha coefficient exceeds the horizontal coefficient, 
and has the opposite sign. If rotation is stronger, the ver- 
tical a effect is strongly reduced. 

For isotropic turbulence, the first-order smoothing ap- 
proximation (FOSA) yields that a can be described by a 
single scalar 



1 



-T{U)-U) 



(29) 



1980D 



where r is the correlation time of the convectio n, an d 
(u) ■ u) is the kinetic helicity (Krause & Radler 
An alternative representation, valid under the assump- 
tion of FOSA, relates a to th e current helicity (Keinigs 
1983; Radler & Seehafer 199C; note the extra minus sign 



in Keinigs' definition of a): 

a^-ri{j-b)/\{B)\\ (30) 

It should be noted that ( ^0|) does not reflect the nonlinear 
feedback of the Lorentz force on the flow, but is a purely 
linear result. If one allows for mild anisotropy in the radial 
direction, alpha assumes the form of a tensor, and a de- 
tailed calculation (Steenbeck & Krause 1969) reveals that, 
for slow rotation (Co <C 1), the diagonal term is given by 

16 



15 



2 2 



3r2-Vln(/7Ui.nis) . 



(31) 
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Fig. 3. Row I : vertical a effect. Row 2: horizontal a effect. Row 3: kinetic helicity. Row 4 ■ gradients of In p and In Urms, 
and their sum (rows 3 and 4 are for a vertical imposed field, but since the field is weak the results in these rows would 
be the same for a horizontal field). Row 5: current helicity for a vertical imposed magnetic field. Left column: run 6 
(Ta = 2000). Middle column: run 7 (Ta = 10 000). Right column: run 10 (Ta = 10 000; no stable layers). In aU cases 
the unstable region is between 2 = and z ^ 1. All calculations are for the solar south pole. Error margins are given 
for rows 1 and 2 as thin lines enclosing the mean; those in rows 3-5 would be comparable or narrower, but are omitted 
in order not to overload the drawings. 



While conditions under which expressions (^9|)-(|3^) hold 
are fulfilled neither in the simulations nor in the Sun 
(mainly because of the requirements for a short correla- 
tion time and only mild anisotropics in the turbulence), 
it is still instructive to compare them with the numerical 
results. 



The kinetic helicity is positive at the top of the convec- 
tive layer (Fig. |^, row 3), as is expected for the southern 
hemisphere based on the direction of the Coriolis force for 
contracting sinki ng pa rcels, or expanding rising parcels. 
Brummell et al. ( 1998 ) find kinetic helicity with the same 
sign for a box on the northern hemisphere, but their coor- 
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dinate system turns out to have a left-handed orientation. 
Near the center of the convective layer, the sign of the 
kinetic helicity reverses if rotation is sufficiently strong, 
because below a certain point sinking parcels expand lat- 
erally while approaching the level z = 1 (i.e., the sta- 
bly stratified overshoot layer or the impenetrable lower 
boundary if an overshoot layer is absent), while rising 
parcels contract while moving away from it. For Co ^ 2, 
no reverse-helicity region seems to exist. Nevertheless, in 
these cases alpha still changes its sign near the bottom of 
the convective layer. For comparison we have calculated 
a case (run 10) with the same Taylor number as in run 
7, but where the stably stratified layers are replaced by 
impenetrable boundaries at z = and z = 1. In this case 
the kinetic helicity near the bottom of the unstable layer 
is more strongly negative, and has a magnitude compara- 
ble to that in the upper part of the unstable layer (Fig. ^, 
row 3, col. 3). This can be attributed to the fact that an 
impenetrable boundary forces a sinking parcel to diverge 
more strongly than does a stably stratified layer. Thus, 
except for the stable region in some cases, the sign of a 
as predicted by Eq. ( p9| ) roughly agrees with that of an, 
while it is opposite to that of ay- The deviating sign of 
ay was explained by Brandenburg et al. ( 1990[ ) . We note 
that the same unconventional sign of ay was also found 
by Ferriere ( |1992| ), and Riidiger & Kitchatinov ( |1993D , al- 
though not for the same values of the Coriolis parameter 
as in our simulations. 

The gradient of In (/OWims) is positive at the top of the 
convective layer and changes sign at a point near the bot- 
tom of the unstable layer in the runs with an overshoot 
layer (Fig. ^row 4, cols. 1,2). Hence the sign of a as pre- 
dicted by (3l|) agrees with that of ajj in the simulations, 
also if rotation is weak. In this respect, V In (purms) agrees 
qualitatively more closely than does the kinetic or the 
current helicity. In the upper region, V In (puims) is pos- 
itive because the density gradient dominates, while near 
the bottom of the unstable domain the influence of the 
convective stability is felt, resulting in a reduction of the 
turbulent velocity sufficient to produce a negative sign. 
With regard to ay, the sign and the dependence on Co 
differ strongly from (|3l|). Also, in the run without the sta- 
bly stratified layer, the gradient of the density dominates 
throughout the box, and no sign change of Vln(pui.nis) 
is observed (Fig. ^, row 4, col. 3). Both components of 
a do exhibit a sign change, though. This inconsistency 
may be a result of the impenetrable boundary condition 
at z = 1 , which causes a transfer of kinetic energy from 
the vertical to the horizontal components, an effect which 
is not accounted for by Eq. (^l|). A detailed comparison 
of the rotational dependence of ay and an with Eq. (^l|), 
and with the analytical results of Riidiger & Kitchatinov 
( |1993D , wiU be presented in Sect. 4.2. 

The current helicity is negative in the upper part of the 
unstable layer for a vertical imposed magnetic field (Fig. 
^, r ow 5) , which confirms earlier results of Brandenburg et 
al. ( 1990 ). A reverse-polarity layer is not observed, except 
in run 10 with the impenetrable boundary at z = 1. Thus, 



the sign of a in the unstable layer as predicted by Eq. (30) 
is the same as that of ay. The negative sign of (j-b) , which 
corresponds to a positive sign in the bulk of the convection 
zone on the northern hemisphere, is at odds with forced 



isotropic turbulence calculations by Brandenburg (2001) 
as well as with observations of the solar surface (Seehafer 
1990, Low 1996) and mean- field calculations (Radler & 
Seehafer 1990| , Riidiger et al. 2001). For a horizontal im- 
posed magnetic field, the current helicity is about one or- 
der of magnitude smaller than in the case of a vertical 
imposed field. It is highly fluctuating, and its mean value 
can have either sign (not shown). The smaller amplitude 
is consistent with the smaller amplitude of the magnetic 
fluctuations (Fig. ||). 

4. The a Effect as a Function of 

In this section we investigate the dependence of the a ef- 
fect on the angular velocity SI. The relevant dimensionless 
parameters are the Taylor number and the Coriolis num- 
ber, as defined by (22) and (23). Averages are calculated 
over time, and separately over the unstable region and the 
stable region, as explained above. 

4.1. Kinetic Helicity and Current IHelicity 

Figure ^ shows the rotational dependence of the kinetic 
helicity and the current helicity (for the vertical imposed 
magnetic field). The kinetic helicity and the current he- 
licity have been multiplied with d/Su^.^-^^ and tj/BqUi-^s, 
respectively, in order to allow a comparison with a/umis, 
as suggested by Eqs. ( p9|) and (|o|) if one sets r « d/urms- 
The kinetic helicity in the unstable layer increases with 
the rotation rate up to Ta « 10"*; its value in the stable 
layer is much smaller. Hence the general trend and the cor- 
rect sign of an (Fig. 5, bottom) are recovered by (^9|) in 
the unstable layer, but not in the stable layer. The current 
helicity for the vertical imposed magnetic field increases 
with rotation up to Ta « 100 — 300, while it decreases for 



0.0 0.04 0.09 0.16 0.30 0.50 



1.34 3.295.38 10.7 




3 



10 



log Ta 

Fig. 4. Solid: Kinetic helicity, {u-u>)d/3u'^^^; dotted: cur- 
rent helicity, r/ {j -b) / B^u^^s, for a vertical imposed mag- 
netic field at the solar south pole. (a,c) average over 
—0.15 < z < 1; (&,c!) average over 1 < z < 1.5. 
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stronger rotation, as do the magnetic fluctuations (Fig. 
It seems that ay in the unstable layer roughly traces the 
current helicity (Figs. |4| and 5). Note however that both 
the normalized kinetic helicity, (a;-it)/a;rms'"rms, and the 
normalized current helicity, (j-b)/jrms^rms, increase mono- 
tonically with Ta (not shown) . This reflects an increasing 
degree of alignment of vorticity with velocity, and of cur- 
rent with magnetic field, respectively. 

The main conclusion to be drawn with regard to the 
validity of (p9|)-( |3l|) is that qualitative agreement exists 
only between an, (|29|), and (^) on the one hand, and be- 
tween ay and (|30| ) in the unstable layer on the other hand. 
For Co ^ 5, such qualitative similarities are no longer ev- 
ident, but in the following sub-section we shall compare 
the behavior of ay and ajj with mean-field results that 
were derived for the case of arbitrary rates of rotation. 

4.2. Rotational quenching of a 

Figure 5 shows ay and an, normalized by itrms, as a func- 
tion of Ta (or Co). For zero rotation the a effect van- 
ishes. For weak rotation ay/urms increases with increas- 
ing rotation rate. It reaches a maximum near Co w 1.3 
(Ta w 2 • 10'^); for stronger rotation it is quenched. The 
decrease of ay/urms occurs in spite of the monotonic in- 
crease of the kinetic helicity, and in spite of the decrease of 
Urms with increasing rotation rate (Fig. The horizontal 
a effect sets in at a higher rotation rate. Its magnitude 
increases monotonically with the rotation rate, and ex- 
ceeds that of the vertical a effect if Co 8. For Ta = 10^, 
Coriolis forces are so strong that they stifle the convec- 
tion (see below), and hence also alpha. Two effects that 
are not shown in detail appear to be responsible for the 
reduction of ay: an increasing number of sign changes of 
ay(z), resulting in cancelations when the volume average 
is calculated, and a decreasing amplitude of the magnetic 
fluctuations (Fig. ||). Both effects are absent in the case of 
an- 

Clearly, the alpha effect is highly anisotropic, and its 
dependence on rotation is more complicated than sug- 

iV). Riidiger 



0.0 0.04 0.09 0.16 0.30 0.50 1.34 3.295.38 10.7 



gested by Eqs. (^ and (|31|). Riidiger & Kichatinov (|993|) 
present analytical expressions for the rotational depen- 
dence of ay and ajj, based on the first-order smoothing 
approximation (FOSA), but for arbitrary rotation rates. 
They employ a linearized equation of motion for the fluc- 
tuating quantities which includes a random driving force 
with prescribed spectral properties to model turbulence, 
and which takes into account a stratification of the den- 
sity and of the turbulent velocity. Due to the stratifica- 
tion as well as rotation, the electromotive force becomes 
anisotropic. In the simplest case of a driving force with 
zero frequency and containing a single wavelength (the 
case designated mixing-length approximation by the au- 
thors), and for a weak magnetic field, a manageable ex- 
pression can be derived. For the present geometry, i.e., a 
local Cartesian grid situated at the south pole, one obtains 



ay 



i/'y V In p ^ 



'i/'yV In Ur 



(32) 




^log Ta 



Co 
1.32 



3.34 5.41 11.6 
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-0.02 
-0.04 
-0.06 



Ruediger & Kichatinov 1993 

«H,slob/ Ur, 




2.0 



2.5 



3.0 



3.5 4.0 



4.5 



log Ta 



Fig. 5. a-coefficients, normalized by Urms, as functions of 
Ta. Here Urms is an average over time and over the unstable 
layer. Top: vertical a effect. Bottom: horizontal a effect. 
The errors are standard deviations divided by the square 
root of the number of turnover times covered by the sim- 
ulation. All simulations are for the solar south pole. The 
dashed curves are analytical results based on Riidiger & 



Kichatinov (1993). In each figure, a suitable uniform scal- 
ing factor was applied to the analytical curves. 



an = — Tu 



V'(^Vlnp + V^Vln^ 



(33) 



These expre ssions correspond to those of Riidiger & 
Kichatinov ( 1993 ) if one sets ay = ctzz, Q^h — o:xx, 
ip^ = Co^-P, and = Co*" (cf. their Eq. 3.11). The 
functions ■f/'y" given by 



1 



1 



23 

T 



(a;2 



1) 



2{x^ + 1) 



6 



■ arctan x 



17- 



X 

{x^ - 1)^ 
a;2 + 1 



■ arctan x 



(34) 



(35) 
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Fig. 6. Strength of the velocity and magnetic field pertur- 



bations as a function of Ta, in units of \/gd, and y/fioPogd, 
respectively. The vertical and horizontal orientation of the 
imposed magnetic field are used for ay and an, respec- 
tively. 




Fig. 7. Functions V'v H' which describe the Coriolis num- 
ber dependence of contributions to ay and an due to den- 
sity and turbulence stratification, respectively. 



where x = Co. Several features relevant for our discussion 
should be noted. As is shown in Fig. ^ the functions V'h" 
are almost equal and, for Co ^ 3, they are linear in Co, so 
that Eq. ( |3^ ) reduces to Eq. apart from a factor of 

order unity. Secondly, ipy", and therefore ay, are subject 
to rotational quenching. This unexpected behavior, which 
is not predicted by Eqs. (|2g|)-(|3l|), is roughly confirmed by 
the simulations, as is shown by Fig. 5. In order to produce 
the dashed curves, representative values, estimated from 
the simulations, were inserted for the gradients V Inp and 
VlnUrms, namely 1.5 and —1.0 in the unstable layer, and 
1.0 and —3.5 in the stable layer (Fig. ^, row 4). For a given 
Taylor number, the same value of Co was used for both the 
unstable layer and the stable layer, i.e. the turnover time, 
T « ^/urms, was assumed to be the same. This is motivated 
by the fact that both the effective thickness, £, and Urms 
for the stable layer are smaller by a factor 3-4 compared to 
their values in the unstable layer. Secondly, the temporal 
variations of the fiow are very similar in both layers (Fig. 



H, bottom). After dividing Eqs. (32) and ( ^ ) by Wrms, a 
factor TUi-ms remains in front which can be assumed to 
be independent of Co. However, it was found that setting 
this factor to rwims = ^ overestimates the amplitude of a. 
For the dashed curves shown in Fig. 5, it was replaced by 
a suitable uniform scaling factor, namely 0.19 for ay and 
0.01 for an- Thus, our simulations produce much smaller 
a coefficients than predicted by ( |3^ ) and (33), especially 
an which is smaller by two orders of magnitude. On the 
other hand, the main point is the agreement with regard 
to the rotational dependence. Except for the sign of ay in 
the unstable layer, the rotational quenching is reproduced, 
including the position of the maxima. For weak rotation, 
ipy « 4?/'y, SO that the contribution of Vlnumis domi- 
nates, while for Co <^ 2 the density gradient dominates; in 
both cases the mean-field formula yields the wrong nega- 
tive sign. Clearly, the sign of ay as predicted by Eq. (32) 
can have both values in principle, and depends in a deli- 
cate way on the relative importance of the two gradients; 
even a slight change of the V' coefficients can generate a 
sign chance. With regard to ajj, we confirm the saturation 
predicted by Eq. ( ^3| ) only in the stably stratified layer, 
but not (yet) in the unstable layer (Fig. 5, bottom). Apart 
from this, the general trend as well as the correct signs are 
reproduced by (^) in both layers. Due to the shape of the 
functions V'h"' ^^^^ feature should be rather robust. 

In order to shed further light on the infiuence of rota- 
tion. Fig. ^ shows snapshots of the vertical velocity for two 
simulations with different rotation rates. Clearly, stronger 
rotation causes the eddies to have a smaller diameter in 
the horizontal plane (i.e. perpendicular to the rotation 
axis). Hence the number of eddies within the simulation 
domain increases with the rotation rate, and this leads to 
a more accurate result for a, as is obvious from the length 
of the error bars in Fig. 5, taking into account the duration 
of the runs (Table |^). In the simpler case of Boussinesq 
convection in a rotating layer, the scale at which the insta- 
bility sets in also decreases with increasing Taylor number 
(Chandrasekhar 1961 ). Although the present simulations 
are done at supercritical Rayleigh numbers whereas the 
analysis of Chandrasekhar applies to the marginally sta- 
ble case, we may for the moment set aside this difference, 
as well as possible effects of compressibility and different 
boundary conditions. In fact, the smallest value of the 
local Rayleigh number in the unstable layer was found 
to be typically lO'', so that the actual supercriticality is 
less than assumed. The present runs, which are charac- 
terized by Pr = Pm = 1 and Ch = 0.018, belong to a 
regime where the onset of instability occurs in the form of 
non-oscillatory convection. Since the value of the critical 
Rayleigh number depends on the horizontal wavenumber, 
k, there is a critical wavenumber, k^, corresponding to a 
preferred eddy size 27r/fcc, for which the critical Rayleigh 
number attains a minimal value, Rac. The eddy sizes ob- 
served in the simulations are in quite good agreement with 
the analytical values of 27r/fcc for convection between two 
free boundary surfaces (Table ||). The results suggest that 
the preferred mode of convection is the same as that per- 
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Table 2. Eddy sizes and preferred length scales for the onset of the convective instability between two free boundaries 



as given by Chandrasekhar ( 1961 ). The eddy sizes were determined visually; their random error is estimated to be 
about 25%. 



run 


1 


2 


3 


4 


5 


6 


7 


8 


9 


Ta 


2 


10 


30 


100 


300 


2000 


10* 


2 ■ 10* 


5-10* 


eddy size 


2.5 


2.5 


2.5 


2.5 


2.0 


1.5 


1.3 


1.0 


0.8 


27r/fcc 


2.8 


2.8 


2.7 


2.4 


2.1 


1.5 


1.1 


1.0 


0.82 




Fig. 8. Snapshot of the vertical velocity (connected lanes: 
downward, isolated patches: upward). Top: run 6, moder- 
ate rotation. Bottom: run 9, strong rotation. The white 
curve denotes the zero level of the vertical velocity. The 
top and bottom surfaces correspond to z = 0.09 and 
z = 1.81, respectively. 

taining to the onset of convection. The critical Rayleigh 
number increases with Ta from 661 for Ta = 2 up to 
2.13 • 10"* for Ta = 10^. This explains why for a fixed value 
of Ra the intensity of convection, as measured by Wrms, 
decreases with increasing rotation rate. Since the actual 
Rayleigh number in the unstable layer has a local mini- 



mum of about 10^, this also explains why no convection 
was observed for Ta = 10^. 



Fluxes 




0.0 0.5 1.0 1.5 



z 



Fluxes 
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z 

Fig. 9. Energy fluxes. Top: run 6; Bottom: run 9. Only the 
three most important contributions are shown, the viscous 
and electromagnetic fluxes being negligible. The total flux 
Ftot does not include the flux associated with the cooling 
term Q, Eq. (||); hence the steep rise in region 1. 

As the scale of the individual eddies decreases, the 
corresponding decrease of the effective Coriolis number, 
251^/ui.ms, may be interpreted as a readjustment of the 
convection in response to increased Coriolis forces, which 
tend to hamper the convective energy transport. Still the 
convective flux of the run with the strongest rotation (run 
9) is reduced by a factor of about two compared to runs 
with weak rotation. Apart from a possible effect of the 
increased Coriolis force, the reduced flux may result from 
the fact that smaller eddies lead to larger temperature gra- 
dients in horizontal planes, and this causes increased hor- 



12 



M. Ossendrijver et al.: a effect: dependence on rotation and magnetic field 



0.08 
0.06 
0.04 
0.02 



0.00 



horizontal 



Bn vertical 



0.0 



0.5 1.0 

Bo/Be, 



1.5 



0.10 
0.08 

0.06 

0.04 

0.02 
0.00 



0.0 



0.5 



1.0 



1.5 



Bo/B 



eq 



Fig. 10. Top: rms velocity perturbations as a function of 
the strength of the imposed magnetic field (normalized to 
the value of Beq for small Bq). Bottom: rms magnetic field 
and final equipartition field strength. 



izontal radiative transport between adjacent eddies which 
reduces the efficiency of convective energy transport in the 
vertical direction (Fig. ||). 

5. Alpha and magnetic field strength 

If the strength of the imposed magnetic field is in- 
creased, effects of the Lorentz force become more notice- 
able. Moreover, these effects depend on the orientation of 
the imposed field, i.e. vertical or horizontal. From the pre- 
viously discussed runs, a representative case (run 7) was 
selected as the starting point for runs with increasingly 
strong imposed magnetic fields of both orientations (Table 
^) . A comprehensive study of the mode structure of mag- 
netoconvection for different imposed magnetic fields is be- 
yond the scope of the present paper, though. Various cases 
for a vertical imposed field, without rotation, are treated 
by Matthews et al. ( |1995D . Wissink et al. (2000) consider 
the breakup of a horizontal magnetic layer into a number 
of flux tubes, thereby including the effect of rotation. 

For an imposed magnetic field in the vertical direction, 
convection is increasingly hampered with increasing field 
strength, as is indicated by a decrease of the rms turbulent 
velocity (Fig. |l^, top) and of the convective energy flux. 
Since the magnetic energy is enhanced by advection of the 
imposed field, while the kinetic energy decreases as a result 
of the Lorentz force, after a while the rms magnetic field 



Table 3. Initial field strength, Bq, and averaging inter- 
val, T. Other parameters as in run 7 (Table |lj). The im- 
posed magnetic field is Bq = B^e^, or BqBx, depending 
on whether ay or ajj is determined. 



run 
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7a 


7b 


7c 
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7e 


So 


0.001 
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T (qv) 
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Fig. 11. Top: av.unst and ay.stab as functions of Bo/Bcq 
for run 7, where Bq = BqEz- Bottom: an.unst and an, stab 
as functions of i3o/^eq for run 7, where Bg — BqE^- Here 
Bcq = 0.052 is the equipartition field strength, averaged 
over the unstable layer, for convection without a magnetic 
field. Analytical quenching results [dash-dotted) are given 
for both signs. All calculations are for the solar south pole. 



can be much larger than the equipartition field strength 
(Fig. bottom). Moreover, if the initial magnetic field 
is too strong, then convection dies out completely. This in 
fact occured during run 7e (Bo/Bcq « 1.7). After 100 time 
units, the volume- averaged energy density of the velocity 
and magnetic-field perturbations had both fallen by about 
2 orders of magnitude. The decrease of ay is more rapid 
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than that of Urms, because advection of the imposed field 
by convection is stifled; this explains why alpha decreases 
also when measured in dynamical units (i.e., when divided 

by Mrms). 

A debate is ongoing in the literature about how and 
whether the magnetic Reynolds number affects the mag- 
netic quenching of the a effect. The dependence of a 
on the magnetic field strength is often schematically 
represented in the form a w q;o[1 + Rm^B'^/B^q]^^. 
While some numerical and analytical results suggest that 
w » (Kraichnan 1979a, 1979b, Brandenburg & Donner 
1997 ), others suggest p « 1 (Catt aneo & Vain shtein 
Vainshtein & Cattaneo 1992, Tao et al 



1991, 



1993 



Cattaneo 1994 Cat taneo & Hughes 1996, Vainshtein 



1998 



Brandenburg |200l|) . Although the present set of simula- 
tions is not well-suited to answer this question because 
of the relatively small Reynolds numbers (Rm « 30), 
the results for both ay and an are consistent with p k, 
1 (Fig. |ll|). Blackman & Field ( |200C| ) argued that the 
quenching with p w 1 observed in simulations with pe- 
riodic boundaries may in fact be a consequence of the 
boundary conditions rather than a dynamic effect, and 
that p = should be possible if a flux of magnetic helicity 
through the boundaries is allowed. Our conditions (|l^) do 
allow for a flux of magnetic helicity through the top and 
bottom boundaries. Thus, the fact that a quenching sets 
in when ^ Rmi?^q suggests that an Rm-dependence 
of the a quenching is possible in the present case; such 
quenching could be dynamic according to the reasoning of 
Blackman and Field. However, this conclusion may change 
if one allows for shear near the surface. This is the case 
between the solar surface and the corona, and so a strong 
helici ty flux is possible in that case (Berger & Ruzmaikin 
2000|) . 

It appears that for boundaries that allow for a flux of 
helicity diverse Rm dependencies of the a quenching are 
possible, depending, e.g., on the character of the forcing 
of the flow. Brandenburg and Dobler (2001) report a case 
where p w 1/2. We also note that even with the strong 
p « 1-quenching of a a large-scale field may build up 
because of an equally strong quenching of the turbulent 
magnetic diffusivity (Brandenburg 2001). 

The magnetic-field dependence of the a coefficients 
is essentially the same for the cases of the vertical and 
horizontal imposed fields, even though convection is ham- 
pered less in the case of a horizontal imposed field (Fig. 
pi] , top). For field strengths comparable to the equiparti- 
tion value, the convective pattern becomes very different. 
If Bq — 0.02 w 0.38Boq, convection exhibits a pattern of 
irregular oblique elongated cells that make an angle with 
the direction of the imposed magnetic field, and a small 
horizontal a effect is still measured. For runs with field 
strengths = 0.04 « 0.77Beq and Bq = 0.06 w 1.15Seq 
respectively, longitudinal rolls are formed that are increas- 
ingly aligned along the x-direction, and the coefficient an 
is practically zero. For still stronger fields [Bq — 0.08 « 
1.5Beq)j convection assumes the form of highly regular 
oblique lanes, and a small a effect is again measured. 



6. Discussion 

The following discussion of our results in the solar con- 
text may be considered as complimentary to that of 
Brandenburg et al. ( 1990| ). First-order smoothing is a jus- 
tified approximation if either the magnetic Reynolds num- 
ber is small, or if the flow is slow in the sense wr/d <C 1. 
In our dimensionless variables, the latter condition would 
mean u-caxsT/d <C 1 but, with the rms velocities and co- 
herence times obtained in our simulations, we rather have 
UrmsT/d 1. Our magnetic Reynolds number (20-35) 
is quite moderate but not small. Hence the simulations 
clearly go beyond first-order smoothing, although we are 
far from the conditions met in the solar convection zone. 
The large fluctuations that occur at large Reynolds num- 
ber are evident in our calculations. In particular we see 
that the components of the a tensor undergo large fluc- 
tuations. Only after averaging over many coherence times 
and over the horizontal coordinates of the box a significant 
depth variation is obtained. With this variation we confirm 
some results obtained in the case of weakly anisotropic tur- 
bulence, and extend them beyond the limits of first-order 
smoothing: in the southern hemisphere the horizontal a 
coefficient (Fig. ||, row 2) is negative (positive) in the up- 
per (lower) part of the box, as is plausible from the ro- 
tational effect on convergent and divergent velocities that 
arise as the flow is forced into the horizontal direction. On 
the other hand we find the opposite sign for the vertical 
a coefficient (Fig. ^ and Fig. |[ row 1). Also of opposite 
sign is the small-scale current helicity. Indeed, both the 
sign and the rotational dependence of ay in the unsta- 
ble layer roughly agree with that predicted by the simple 
FOSA-expression which relates a to the current helicity 

&■ 

For cases with convective overshooting in the lower 
part of the box the components of a reverse their sign 
near the transition to the stable stratification. For an aQ 
dynamo the horizontal component an is the essential in- 
gredient; possibly the sign change of an has interesting 
implications for a dynamo that operates as an interface 
wave at the base of the convection zone. In the layer of 
convective overshooting the sign of an is appropriate for 
the equatorward migration of the mean field, if the posi- 
tive shear dQ/dr, as inferred from hclioseismology for the 
depth range around 0.7ro and the lower heliographic lat- 
itudes (Kosovichev et al. 1997), is used as the other es- 
sential ingredient. These conclusions are drawn from sim- 
ulations with an imposed magnetic field Bq ~ 10^^, and 
may still be valid for Bq — 10~^, as the results of run 7a, 
shown in Fig. |ri|, suggest. The latter field strength corre- 
sponds to w 10 T, which is the strength of the field in the 
stable layer below the convection zone, as deduced from 



the properties of rising flux tubes (Caligari et al. 1995). 
For stronger fields we obtain a marked quenching of a con- 
sistent with a factor [1 -I- Km {Bq / Bcq)^]~^ , although we 
have checked this only for Rm « 30. The simulations have 
also revealed a striking difference between the convective 
patterns for a vertical and a horizontal orientation of the 
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imposed magnetic field, with rolls being formed for strong 
horizontal fields. 

In view of a strong toroidal field at the base of the 
convection zone a magnetic Rayleigh- Taylor instability 
(Brandenburg & Schmitt 199^) or an instability of toroidal 
flux tubes (Ferriz-Mas et al. 1994) should be employed to 
obtain an a effect. Nevertheless, since the sign of an fol- 
lows plausible arguments, we may hope to obtain the same 
sign for such a dynamically determined a (Brandenburg 
& Schmitt |1998D . 

As far as the dependence of a on rotation is concerned, 
the simulation may yield more reliable results than for the 
other parameters. In the solar convection zone the Coriolis 
number is small near the surface, and increases to values 
around 1 or somewhat larger near its base. This is the 
regime shown in Fig. 5. The rotational quenching of ay 
and the saturation of ajj (Fig- 5, top) may therefore just 
begin in the depth region where the dynamo operates, but 
in solar-type stars which are younger and therefore rotat- 
ing faster, these effects may be important. We have ob- 
served some agreement between the numerical results for a 
and the predictions of Riidiger & Kichatinov (1993), as far 
as the dependence on rotation and on depth is concerned. 
In both the unstable and the stable layer, the approximate 
Coriolis number dependence and the correct sign for an 
are reproduced. The main discrepancies concern the am- 
plitude of a (the scaling factors of the theoretical curves 
in Fig. 5) and the sign of ay in the unstable layer, as well 
as the sign of the small-scale current helicity. 

Even within the limited range accessible to numeri- 
cal simulation, the parameter space is enormous. We have 
not yet considered the latitude dependence of a (although 
some earlier work of Brandenburg ( 1994 ) indicates that 
\a\ is maximum at about ±60° of latitude), and we have 
not yet explored the dependencies on the Reynolds and 
Rayleigh numbers. This must be done by further simula- 
tions. 
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